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ABSTRACT 



A new more rigorous and accurate method for treating neutrino oscil- 
lations in the context of the MSW effect in a medium is proposed. This 
leads to a new type of resonance condition which for small mixing angles 
puts rather stringent conditions on E y /8m 2 . The implications on the solar 
neutrino problem are discussed. 
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One of the most important questions of the physics beyond the standard 
model is the problem of neutrino masses. Furthermore if the neutrinos are 
massive, the neutral leptons produced in weak interactions are not station- 
ary. They are linear combinations of the neutrino mass eigenstates (neutrino 
mixing). Even though at present there is no theory which can predict the 
mass and mixing of the neutrinos, most Grand Unified Models predict small 
masses and mixing. If the neutrinos are almost degenerate, neutrino oscilla- 
tion experiments o @ are the best candidates to measure small 6m 2 (from 
leV 2 down to 10 _10 eV 2 ). Furthermore, neutrino oscillations may explain 
the solar neutrino problem ® S & S, i.e. the apparent reduction of the 
flux at earth compared to that predicted by the standard solar mo del I 
(SSM). The mechanism of neutrino oscillations however, is not effective if 
the neutrino mixing is small. It has been observed, though, that under the 
conditions of high density encountered in the sun's interior the oscillation 
can be enhanced due to the MSW effect In other words small mixing 

angles can be converted into large effective mixing angles due to the resonant 
scattering of v e neutrinos by electrons. Calculations of neutrino oscillations 
involving the MSW effect have hitherto involved the following steps: 

1) Convert the evolution equation from flavor space into neutrino mass 
eigenstate basis by locally diagonalizing the space dependent HamiltonianS' 

2) Ignore the transitions between the mass eigenstates (adiabatic approx- 
imation) or treat such transitions in perturbation treatment. In the last case 
one assumes that &3. 

6E v /2 



d'dm/dx 



« 1 (1) 



where 5E U is the difference in neutrino energies and the variation of the 
mixing angle in the medium with distance. In particular the above equation 
must be true at the resonance point. 

In the present paper we will provide an exact solution which does not go 
through the local neutrino mass eigenstates. Our method is quite simple and 
offers itself to a simple interpretation. 

For illustration purposes we will exhibit our method in the case of two 
generations but it can easily be extended to any number of generations. 



Following standard procedure tUJ we can write the neutrino state at any 
time as 

\u(t) >= a e (t)\u e > +a a (t)\u a > 
where \v e > is the electron neutrino and v a any other flavor (e.g. v^). The 
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amplitudes a e (t) and a a (t) satisfy the evolution equation 

= H(t) 



d_ 

% dt 



(In 



(2) 



Since x = ct, the above equation can be written in term of x. In the presence 
of matter it can be shown that TC can be cast in the form o 



H 



E u - f cos2$ + ^ 
jsin2 , d 



E v + f cos2d 



where 



(3) 



(4) 



5m? (5m 2 /leV 2 ) 

$ is the usual (vacuum) mixing angle and £o(x) takes into account the fact 
that the charged current interaction between z/ e and electrons causes a shift 
in the electron neutrino energy. £q(x) takes the form 



£ (x) 



Air 



2V2G FPe 



x 



(5) 



where Gf = 1.16636 x 10~ 5 GeV~ 2 is the Fermi coupling constant, and p e (x) 
in the number of electrons per unit volume at distance x from the sun's center 
which is assumed to be spherically symmetric. Equation @) can easily be 
integrated to yield 



exp(— i A) 



(6) 



where 



a n 



is the initial solution of the differential equation and the ma- 



trix A is given by 
A 

where 



TT(X 



Xq) 



2^-cos2d sin2d 
sin2{} cos2d 



x dx' 



(7) 



(8) 



(x - X ) Jxo £q(x') 

The diagonal matrix with elements E v {x — xq) has been omitted as irrelevant. 

Using the well known fact that the minimum polynomial associated with 
the matrix A is of degree one, which permits exp(— lA) to be written as a 
linear combination of the identity matrix and A plus the fact that 



\v(xq) >= \v e > 



2 



it is straightforward to show that 



C2 



2 (X — Xq) 



V(u e - u e ) = 1 - -£sin 2 20sin 2 ir K " (9) 
where 

£m = £ m (x,X ) (10) 



£ 

(i - cos2$) 2 + sm 2 2tf 

£ m can be interpreted as the oscillation length. In the presence of matter, 
however, £ m is a function not only of the combination (E u /8m 2 ), but also 
of x (detector point) and xo (source point). One can also talk about an 
effective mixing angle d m defined by 

sm2i? m (j,i ) - ^ £ — —sin2d (11) 

From eq.flTTD it can easily be seen that the maximum oscillation probability 
occurs when 

£ = cos2d , ("resonance" condition) (12) 

In this case 

sin2$ m (x, xo) = 1, £ m (x,xo) = £/sin2'd (13) 
while eq.(^|) becomes 

P (,.^,J = l-*,. { ,^} I") 



The above condition ( |T2] ) is reminiscent but not similar to the usual reso- 
nance condition (see eqs (|25|) and (|26|) below) which occurs at some appro- 
priate point Xr. In our exact treatment the "resonance" condition (eq. (|T2"D) 
involves both the initial (xq) and the final (x) positions. It is affected by 
the cumulative effect of the density of the medium and not by its value at 
some appropriate point xr. It is a sort of "global resonance" condition. For 
numerical calculations we will cast £ in the form 

£ = f X ^dx' (15) 

(x-Xq) Jxq po 

where 

_ 2y/2G F E vPo 
to ~ 6n^ (16) 
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with po the sun's density at some suitable point (e.g. at its center). The 
gross features of the resonance are affected by to, while its details depend on 
the specific from of p e (x)/po. The value of to which is consistent with the 
present solar neutrino data, lies in the range 1 



0.5 < t < 50 (17) 
This range maybe enlarged however (~ 0.2 — 60), if also uncertainties on the 



density p are also includedl!3 . The resonance condition imposes constraints 
on the parameters Xo,5m 2 and E v (x is assumed to be fixed, i.e. x = sun - 
earth distance). The half maximum width is given by 

T = 2sin2-d (18) 

which for small mixing angle puts stringent constraints on t or equivalently 
on the allowed values of E y /5m 2 . 

Before proceeding further with the discussion of our results we will com- 
pare the above new formulas with those which have been obtained with the 
traditional approach. By diagonalizing the matrix of eq. @ for each value 
of x we obtain the eigenvaluesE H3 

"^^i+uD (19 » 

and the eigenvectors 

y= cos$ m (x)\v e > sin-d m (x)\v a y (20) 

\vh{x) y= siwd m (x)\u e y +cos , d m (x)\v a y (21) 

with 

tW^/H^-^.^jT" (22) 

£ £ 

sin2$ m (x) = sin2d/\l - 2cos2$-— + (^-) 2 } 1/2 (23) 

£q{x) £ {x) 

The resonance in this approach occurs when the diagonal elements of 7i of 
eq.(|3|) are equal 0, i.e. at a point xr in the medium such that 

£ (x R )cos2d = £. (24) 

At the resonance one finds 

sin2 r d m {xRj = 1 (25) 

£ m (x R ) = £/sin2$ (26) 
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Then by writing 

\v(x) >= a L {x)\v L {x) > +a H \v H (x) > 
we obtain the following evolution equation 



% d_ I a L {x) 
dt \ a H {x) 



B{x) 



a L (x) 
a H {x) 



where 



B 



ax 



(27) 



(28) 



(29) 



Following a procedure analogous to that of our new method outlined above, 
the oscillation probability takes the form 



= -{1 + cos2i) m {x)cos2i) m {x ) — — 

2 jj 2 + Y 

+ sin2$ m (x)sin2$ m (xo)cos(2^J [3 2 + 7 2 ) 

+ sin(2& m (x) - 2$ m (x )) m ^= stn{2J(3 2 + 7 2 )} (30) 

VP +7 



where 
and 



7 = i? m (a;) - fi m (xo) 
r dx 

P = 71 



If the sun's density is discontinuous at its surface, we get 



P = 7T 



1 



(31) 
(32) 

(33) 



where x s = i? is the sun's radius. We notice that the quantities $ m (x) (and 
7) and £ m (x) are not independent quantities but they are given by eqs(|25| 
and |26|). In the special case 7 <C (3 (adiabatic approximation) we obtain 

V[y e — > u e ) = -{1 + cos{2 , d m (x) — 2^ m (x ))} — sin2$ m (x)sin2 , d m (x)sin 2 j3 

2 (34) 
which coincides with the old expression© . Furthermore in the absence of 
matter effects i.e. when $ m (x) = $ m (xo) = $ we obtain once again the well 
known formula 



V{y e — > v e ) = 1 — sin 2 2$sin 2 { 



X — Xq 



7T- 



} 



(35) 
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On the other hand if i? m (x) = $ and sin2$ m (xo) = sin2d(xn) = 1 (i.e. the 
resonance condition occurs at Xq) one obtains for small values @ of $ 

V{v e -> v e ) ^ - - sin2$sin 2 f3 ~ - (36) 

The comparison of our treatment with the non-adiabatic treatment of the 



earlier perturbative calcultations is not obvious. Level crossings E 
etc do not enter in our non-perturbative treatment which avoids the inter- 
mediate step of the eigenstates of eqs.( [20"l - [2~ll ). We only notice that in the 
limit 7>/5 our expression (^) becomes 

V{y e ^u e )~l (37) 



which means no oscillations in this limit(see also eq(91) of ref.||12||). 

Returning back to our new method we repeat that in the case of small mix- 



ing angles for a resonance to occur, eqs (p!2| , |l9|) impose a stringent condition 
on the properties of the neutrinos. Since in the solar neutrino experiments 
the sun-earth distance x is fixed we expect that the resonance will occur for 
special values of xq. To test this we will consider a reasonable model for the 
solar electron density employed by Lim and Marciano GJj, i.e. we take 

(38) 



Pe(x) = f l-ajf, 0<| x |< hR e 

Po I (1 — aki) exp{c{k\ — jr 1 )} \ x \> kiR Q 



with ki = 0.2, a = 10/3, and c = 100/9. The value of po, absorbed in the 
definition of t (see ITST) was chosen 



Po 



6 x 10 25 cm~ 3 (39) 



With the above parameters we have perfomed calculations of V[y e — > u e ) 
as given by eq.(|9|) near the resonance (|T2"D for various values of xq. The 
obtained results for sin 2 2d = 2m e /m fJi » 1.1 x 10 -2 and Xq = 0, Xq = 0.1 R Q 
and Xq = 0.2R Q are shown in figs (1-3). The corresponding values of E u /Sm 2 
in units of MeV — eV~ 2 are also indicated. It is clear that the oscillation 
probability becomes negligible for regions away from the resonance. Our 
results are similar to those obtained for x — in the context of the earlier 
treatment by Rosen and Gelb except that: 

i ) The value of E v j 5m 2 at resonance is in our approach about a factor of 
2 smaller. This, however, can be accounted for by the fact that they use a 
different electron density function. 

ii) The width of our plots of oscillation probability as a function of 
E u /Sm 2 is quite a bit narrower. 
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We are currently analysing the solar neutrino data employing our new 
formalism for a number of solar density profiles. Detailed results will be 
published elsewhere. At present we note that for suitable values of E v /8m 2 
the resonance condition for a small mixing angle is satisfied in a small region 
in the x - space. One must appropriately integrate over x 0. If one does 



this, one expects that even though some neutrinos may arrive at earth with 
maximal oscillation probability as given by eqs ([L4]) and ([13]) the fraction of 
such neutrinos will be small. A rough estimate can be given as follows 

N v (res.) JP(z/ e -> is e )p(x )dx 
N„(tot.) Jp(x )dx 

which in the interesting case of small mixing, i.e. xsin2-d <C I one finds 

NJres.) _ x 2 



N v {tat.) 



-sin 2 2$ (41) 



If this is borne out by more detailed calculations in which integration over 



the neutrino energy is also preformed 113 it will imply that in fact the 
medium may have a negligible effect on neutrino oscillations. 

In conclusion we have presented a novel way of treating neutrino os- 
cillations in a medium which is both simpler and more accurate than the 
traditional perturbative approach. Our method can be easily extended in 
the case of the three generations. Our results given above are not expected 
to be drastically altered by such an extension. 
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Figure Captions 



Figure 1: Oscillation probability at Earth for a neutrino created at the 
center of the sun (x = 0), cl) clS cl function of the neutrino energy E v devided 
by 5m 2 = ml a — m1 e and, b) as a function of the distance from the sun. 



Figure 2: same as in figure 1 but for x = 0.1i? o 



Figure 3: same as in figure 1 but for x = 0.2R Q 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404350vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9404350vl 
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